Entanglement and the nonlinear elastic behavior of forests of coiled carbon nanotubes 
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Helical or coiled nanostructures have been object of intense experimental and theoretical studies 
due to their special electronic and mechanical properties. Recently, it was experimentally reported 
that the dynamical response of foamlike forest of coiled carbon nanotubes under mechanical impact 
exhibits a nonlinear, non-Hertzian behavior, with no trace of plastic deformation. The physical 
origin of this unusual behavior is not yet fully understood. In this work, based on analytical models, 
we show that the entanglement among neighboring coils in the superior part of the forest surface 
must be taken into account for a full description of the strongly nonlinear behavior of the impact 
response of a drop-ball onto a forest of coiled carbon nanotubes. 



The study of nanostructures, in special carbon nano- 
tubes (CNTs) and nanowires have been object of intense 
experimental and theoretical investigations due to their 
large range of possible applications and new physical phe- 
nomena [TJ |2] . Among these nanostructures helical and 
coiled structures have a special place due to their differen- 
tiated mechanical behavior [5] . Coiled carbon nanotubes 
(CCNTs) were first predicted to exist in the early 1990s 
by Dunlap |4] and Ihara et al. [3[6j|7] and experimentally 
observed in 1994 by Zhang et al. [S]. CCNTs have been 
receiving increasing interest because of their additional 
capability to serve as nanoscale mechanical springs [5], 
electrical inductors [10] . and for their potential applica- 
tions in composites [TT] . 

Recently, the dynamical response of a foamlike forest of 
CCNTs (Fig. 1(a)) under impact of a drop-ball has been 
reported |12j . The experiment consisted of producing 
arrays of bundles of CCNTs [T31 [2] , let a stainless steel 
bead falls down on the forest of CCNTs, and measure 
the dynamic force at the wall below the forest during the 
stages of penetration and restitution. The analysis of the 
forest's morphology after impact has shown no trace of 
plastic deformation and a full recovery of the foamlikc 
layer of CCNTs under various impact velocities. 

The contact force exhibits a strongly nonlinear depen- 
dence on displacement and appears fundamentally differ- 
ent from the response of a forest of CNTs [23 EE1 EZ1 EH] ■ 
The obtained results in [T^] have been compared to the 
Hertz elastic model [T3] of a solid sphere and a planar sur- 
face in contact where, if F is the force of contact and S 
is the displacement, F ~ 6 [19]. The results reported 
in Ref. [T^] showed a nonlinear response of the CCNT 
forests having a force-displacement relation of F ~ S 2 2 , 
different than the Hertzian case. It has been hypothe- 
sized that sideways interactions of the compressed CC- 
NTs can be associated with the strong nonlinear behavior 
in CCNT forests, but the physical mechanisms behind it 
are still unclear. 

In this work we propose a model for the impact res- 
ponse of a forest of CCNTs. The model takes into ac- 




FIG. 1: (Color online) (a) Scanning electronic microscopy pic- 
ture of a CCNT forest [T2]. (b) Schematic diagram of a CCNT 
forest deformed by a drop-ball. H is the forest thickness, R is 
the ball radius, and 8 is the total depth displacement of the 
ball into the forest. A is the axial deformation of the spring at 
a distance r from the center of the contact between the ball 
and the forest, (c) Parameters of a helical spring: a is the 
coil radius, P is the spring pitch and D spring diameter, (d) 
Schematic representation of the entanglement between two 
adjacent coil bundles forming inactive, or "solid", and active 
turns. 



count: i) the individual elastic contribution from each 
CCNT in contact with the drop-ball; ii) the geometry of 
the surface of contact between the forest and the drop- 
ball; and iii) sideways interactions of the compressed CC- 
NTs through an entanglement process. As discussed be- 
low we show that this model can describe the strongly 
nonlinear behavior observed in the recent experiments 
with CCNT forests Q2]. 

In the model the CCNT forest of thickness H is consid- 
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ered as a set of identical springs that individually interact 
with the drop-ball and that may laterally interact with 
other springs. The drop-ball is considered to have a rigid 
spherical surface of radius R. Figure [l|b) illustrates the 
penetration of the ball into the CCNT forest and the 
geometric features of the model. The total depth dis- 
placement of the ball into the CCNTs forest is 5. Each 
bundle of CCNTs is modeled as a helical spring of pitch 
P, and with a Hooke's constant k = (G D 4 ) / (64a 3 N) pO] . 
where G is the shear modulus, N is the number of active 
turns, a is the coil radius, and D is the bundle diameter 
(Fig. 1 (c)). Typical values from experiment are H = 
100 fim, R = 1000 urn, 2a = 0.45 fim, D = 0.1 /im, S = 
3 jum, and P = 0.9 /im [12j [14]. We also considered the 
situation where the bundles of CCNTs are very close to 
each other. From Ref. [T5] the total density of CCNTs in 
the forest is ^100//im 2 . For a rough estimative, we ob- 
tained a separation between the centers of adjacent coils 
of ~0.1 yum, assuming an ordered squared grid of CCNTs 
distributed on the surface with the value of the total den- 
sity of CCNTs. This separation value is of the order of 
the estimated bundle diameter (~0.1 /im). Therefore, 
it is reasonable to assume the existence of a certain de- 
gree of entanglement between some parts of the CCNTs 
throughout the forest which would contribute to sideway 
interactions. Furthermore, it is probable that changes 
in the entanglement occur during the dynamic contact 
between the ball and the forest surface. This mecha- 
nism is illustrated in Fig. [ljd), where it is shown that 
the impact of the ball on the forest top surface would 
cause a bending of the tips of each CCNT, leading to 
contacts among the superior turns of the adjacent coils. 
The entanglement mechanism can be very complicated 
and a long-range process. In the present model, the ef- 
fect of this mechanism in the forest response is translated 
into the reduction of the number of active turns of each 
CCNT, i.e., 



forest surface proportional to the ball velocity: 



N = N T 



N s , 



(1) 



where Nt is the total number of turns of each CCNT 
(Nt — H/P = 110), and Ns is the number of turns that 
become inactive due to the entanglement process (Fig. 
1(d)). The inactive turns form a "solid" phase that do 
not contribute to the elastic response of the CCNT. This 
approach has been used by Rodrigues et al. [3T] to obtain 
a nonlinear relation between the force and displacement 
of a conic spring. Based on the small displacement (~ 
3 /im) of the top forest surface compared to the forest 
thickness (~ 100 /im), we also assume a short-range ef- 
fect of the entanglement, i.e., only nearest neighboring 
coils interact. Initial entanglement prior to the impact 
can be incorporated in the initial number of active turns. 
The short-range interaction is included into the model 
by making the rate of increase of the number of inactive 
turns during the contact between the drop-ball and the 



dN s 
dt 



TjV, 



(2) 



where 77 is a measure of how many turns become inactive 
per unit length of the displacement. Despite the com- 
plexity of the interaction that forms the entanglement, r\ 
is considered here constant and determined from experi- 
mental data a 'posteriori. 

Since vdt — dz, where z is the direction perpendicular 
to the forest surface (that is the direction of the move- 
ment of the ball) , Eq. ^ can be easily integrated to give 
the following expression for the number of active turns: 



N(A) = N T - 77A, 



(3) 



where A is the axial deformation of each CCNT. Figure 
[TJb) shows how to determine the value of A as a func- 
tion of 5 and the distance r from the center of the contact 
between the ball and the CCNTs forest (hereafter called 
"center of contact" for short). From geometrical consid- 
erations we have 



A(r) = \[W 



(R-6). 



(4) 



The total force F consists of a summation of the forces 
of each spring that interacts with the drop-ball: 



M 

F = k^ 



1 - ( v /N T )A( n ) 



(•5) 



where i represents the set of rii springs that are at the 
same distance Ti from the center of contact. At the edge 
of the contact area, M is such that A(tm) = 0, i.e., 



M = 



V2RS - 5 2 
2a ' 



(6) 



We assume a circular symmetry of the projection of the 
surface of contact between the ball and the CCNTs forest, 
on the plane of the forest. Therefore, rii is given by 



?io = 1 and rii = 



2jtt\ 
~2r~ 



= 2m, i= 1,2,3, ...M, (7) 



where we used that r, = 2ai. The index i = is for the 
first coil hit by the ball. Equation ^ is, of course, an 
approximation for the situation where the concentration 
of the CCNTs are such that they are beside each other. 
Equation Q will be a good estimative of the number of 
coils at distance r$ of the center of contact if the pene- 
tration of the ball is enough to ensure that many coils 
are hit by the drop-ball. That is the case considered in 
the experiment, where a total circular contact area has a 
radius of ~ 77 /im and the CCNT bundle radius is ~ 0.2 
/i,m [TJ]. For smaller depths Eq. ^ should be corrected 
for a better result. 
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Since S/R -C 1 it is reasonable to use the follow- 
ing approximations: Afr) = 6 — (2i 2 r 2 )/R and M = 
V2R5/(2a). Thus Eq. M can be rewritten as 



F ^ 



kS 



M 



1 - ( V /N T )S 



i(SR-2i 2 a 2 ) 



1 [1~(t 1 /N t )S]R + 2t 1 



(8) 

A simpler expression for F cannot be derived from 
Eq. However, by neglecting the term 2r]i 2 a 2 from the 
denominator, the following expression can be obtained 



F : 



ttRS 2 



4a 2 [l-(r]/NT)S] ' 
where terms proportional to 5/R were also neglected 



(9) 
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FIG. 2: Force-displacement curves obtained from Eq. |8| 
when the term 2rji 2 a 2 is present (solid lines, no approxima- 
tion) or not (dashed lines) for (a) r\ =l/^m, (b) r\ =5//jm, 
and (c) r\ =10//im. The following parameters where used: 
a = 0.225 fim, R = 1000 /Ltm, N T = 110, and k =1 N/m. 



The result of these approximations to the final behav- 
ior of F is shown in Fig. [2] using typical experimental 
parameters. We can see that this approximation is phys- 
ically sound for r/ < 10//im. We will see that this condi- 
tion is satisfied in the experiments [12]. 

Using Eq. ^ to fit the experimental data [T2] we ob- 
tain k = 1.988 N/m and 77 = 7.066//im. The comparison 
with experiment is presented in Figure [3] 

The importance of the entanglement for the nonlinear 
behavior of the forest of CCNTs can be tested by recal- 
culating the total force with r\ — in Eq. ([3]), yielding 
to 



F= k 



nRS 2 



(10) 



Eq. (10 1 shows that the force F is proportional to 8 2 . 
Even being close to the relation found in [T5], this re- 
sult still not capture the full nonlinear behavior of the 
impact response of a forest of CCNTs. However, when 



the entanglement is turned on again, our model, with the 
fitted parameters k and 77, recovers the previous experi- 
mental fitting, F = AS m , A = 0.031 and m = 2.2 [12], 
what shows that the entanglement formed by lateral de- 
formations of the CCNTs is necessary to explain the full 
strongly nonlinear behavior of the impact response of a 
drop-ball onto a forest of CCNTs. 

In a loading experiment, Cheng et al. have measured 
the spring constant value of a single amorphous carbon 
nanocoil in a low-strain regime (nanocoil elongation <3 
/Ltm) as being 0.12 N/m [5]. The number of turns of that 
nanocoil is about 10 which leads to ~1.2 N/m for the 
value of the spring constant of a single turn. Using the 
derived k value of the present model from experimental 
data [12] we can estimate the spring constant of a single 
turn of an individual CCNT. According to Ref. [T^] the 
bundle is formed by ~25 nanocoils. Assuming that the 
bundle is formed by a parallel association of CCNTs, our 
estimative for the spring constant of a single CCNT is 
k s = k/25 = 0.08 N/m. The spring constant of a single 
turn of the CCNT in a compression experiment is, then, 
k s x Nt — 8.7 N/m, which is the same order of magni- 
tude of the value for a single amorphous carbon nanocoil 
measured in a tension experiment. 

As previously mentioned, rj is the number of turns that 
become inactive per unit length of displacement. Sup- 
posing that the entanglement is formed by inactive turns 
that form a "solid " phase in the top of the CCNT forest, 
we can estimate the thickness, t, of the entanglement at 
a given distance r from the center of impact as: 

t(r) = r)DA(r). (11) 

where A(r) is given by Eq. Q. The inset graph of Fig. 
[3] displays t as a function of r. For the center of impact, 
the thickness is 2.1 /im and the number of inactive turns 
of the coil at the center of impact is ~21, about 20% of 
the total number of turns of the total forest thickness. 
These results are compatible with what we expect for an 
entanglement formed by the superior parts of the CCNTs 
due to the impact of a drop-ball on the forest. 

If, instead of a ball, we have an approximated per- 
fect cube or parallelepiped with a finite contact area A 
with the forest, falling down on the forest of CCNTs, 
its impact response can be estimated using Eq. (5). In 
this case, all nanocoils feel the same axial deformation 
<5, and the sum in Eq. (5) can easily performed to 
give F ~ n c S/[l — (rj/N^S], where n c is the number 
of nanocoils in contact with the face of the cube or par- 
allelepiped, and is simply given by n c ~ A/(ira 2 ). The 
response force of the forest, then, will be approximately 
linear with S, for small 6 values, and start growing non- 
linearly for larger 5 ones, as consequence of the entangle- 
ment. 

It should be stressed the limit of validity of our model 
(represented by Eq. ([2])) takes into account the experi- 
mental condition 1121 of small forest deformations. One 
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of the predictions of the model that can be experimen- 
tally tested is the thickness value of the entanglement of 
the top forest surface. We hope the present work will 
stimulate further experiments to test the validity of the 
present model. 




5(|xm) 



FIG. 3: (Color online) Behavior of the force as function of the 
displacemet of the drop-ball during contact with the CCNT 
forest. The model curve was obtained using Eq. |9| with k = 
1.988 N/m and r\ = 7.066//im. The inset graph shows the 
behavior of the thickness of the entanglement as function of 
distance from the center of impact (r = 0). 



In conclusion, we have derived an analytical model for 
the nonlinear behavior of the impact response of a for- 
est of CCNTs including geometrical and physical aspects 
during the forest compression. We showed that the non- 
linear behavior is fully described when the entanglement 
of the coiled carbon nanotubes in the superior part of the 
forest surface is incorporated into the model. This en- 
tanglement among neighbors is due to the bending of the 
coil tips produced by the ball impact. Under the exper- 
imental conditions of small deformations [T^] the model 
predicts an entanglement thickness of ~ 2 /im at the 
maximum forest compression. The model results point 
out to the importance of the coil entanglements for the 
elastic behavior of such systems. The present model is 
able to provide, by matching experimental values, esti- 
mates of the spring constant of a single CCNT and the 
level of entanglement between CCNTs. These aspects 
can play an essential role in the future design of micro- 
electro-mechanical systems devices, new shock protecting 
layers, and composites for microelectronic packaging and 
vibration mitigating materials where CCNT structural 
entanglement could be present. 
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